Systematic study of symmetry energy within the SMM picture of multifragmentation 
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A systematic study on the effect of secondary decay on the symmetry energy coefficient extracted 
by isoscaling and the recently proposed isobaric yield ratio methods within the Statistical Multifrag- 
mentation Model is performed. The correlations between the input symmetry energy coefficients 
and the calculated ones from both primary and secondary fragment yields are analysed. Results for 
secondary fragments show that the best estimation of the input symmetry energy coefficient within 
SMM is obtained by the isoscaling method, using the yields of light fragments. A comparison to 
experimental results is also presented. 

PACS numbers: 21.65Ef, 24.10Pa, 25.70Mn 
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Theoretical predictions [l[ suggest that information on 
the symmetry energy term of the nuclear equation of 
state can be extracted from the isotopic distributions of 
primary fragments produced in multifragmentation reac- 
tions. However, quantitative information is difficult to 
extract as most fragments, produced in excited states 
@j II] i decay to lighter stable isotopes on a typical time 



scale of ~ 10 s [4|, before being detected. These latter 
fragments are commonly referred to as secondary frag- 
ments. Previous work has evaluated the excitation ener- 
gies of primary fragments [3, E| and indicated that sec- 
ondary decay may distort the signatures of the symmetr 
energy contained in primary fragment observables 0. 
It is therefore important to study model predictions for 
observables that can be calculated for both primary and 
secondary fragments. 

A systematic study on the effect of the secondary de- 
cay as predicted by the Statistical Multifragmentation 
Model (SMM) Q is presented in this report, with partic- 
ular emphasis on the comparison between the recently 
proposed isobaric yield ratio method Q and the well 
known isoscaling [y, I9l-[l3j. To ascertain the degree of 
confidence that can be obtained in experimental results, 
the correlation between input values and the quantities 
that can be extracted from secondary fragments was es- 
tablished. A comparison of SMM predictions to exper- 
imental data measured in 78, 86 Kr+ 58 ' 64 Ni reactions at 
35 MeV/nucleon with the NIMROD-ISiS array Q is also 
presented. 
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Isoscaling parameters deduced from isotopic yields 
measured in two similar reactions with different isotor. 
composition are commonly used observables [TTI . 
Ii~7j to access the symmetry energy in heavy-ion colli- 
sions. The statistical interpretation of isoscaling links 
the isoscaling parameter a to the symmet ry energ y coef- 
ficient, C sym , of the equation of state pl1.lla.fl8l]: 
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where T is the 
sources, A = 



temperature of the two fragmenting 
■)j — (3)0 an d the (Z/A)i values cor- 



respond to the proton fraction of the n-poor [i = 1) and 
n-rich (i = 2) sources, respectively [III, EH- Other def- 
initions of the quantity A have been recently suggested 
fl7l - f]~9l . [20I ] . which take into account the fragment isotopic 
composition rather than the source composition. How- 
ever, the debate on the proper choice of A is still open. 
In this work, we will restrict ourselves to the A definition 
suggested in fill . HEl ] , which is the one first used for the 
statistical interpretation of the isoscaling parameters. 

The isobaric yield ratio method, recently proposed in 
Ref. [|[, allows one to extract C sym from the yield ratio 
of two pairs of isobars, A, produced by the same reaction 
systems: 
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where i?(3, 1,^4) and R(l, —1, A) arc the ratio of the 
yields of isobars A, with N — Z = 3, 1 and 1, —1, re- 
spectively. The quantity 6(3, 1,A) is the difference in the 
mixing entropies of isobars A, with N — Z = 3, 1 and can 
be neglected 

In this work the SMM model 0] was used to simulate 
the statistical fragmentation of sources corresponding to 
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the Kr projectile in the reaction systems 78 Kr+ 58 Ni and 
86 Kr+ 64 Ni [2l| that were experimentally measured. This 
will allow the comparison of simulated and experimental 
data. The version of SMM used for the calculation is 
the one described in Ref. [22|, I23I . The break-up den- 
sity was chosen to be p = po/Q [101 1. The input sym- 
metry energy coefficient C™ Tn was varied between and 
the standard value of the symmetry energy at normal 
density, 25 MeV. Source excitation energies (E*) of 3, 5 
and 7 MeV/nucleon were chosen, allowing comparison to 
previous experimental data [21j . To achieve a statistical 
uncertainty on the yield of every analyzed fragment bet- 
ter than 10%, 10 6 events where generated for each case 
analyzed in this paper. Isotopes for which fewer than 100 
counts were obtained in the calculation were not used in 
the analysis. Isotopes with Z > 17 were not included 
in the analysis, since Z = 17 corresponds to the highest 
fragment charge detected with a good isotopic resolution 
by the NIMROD apparatus [Hj], used in the experimen- 
tal campaign. Final fragment yields were obtained by 
using the standard deexcitation procedure implemented 
in SMM. All ground and nucleon-stable excited states 
of light fragments, as well as in-medium modification of 
fragment properties at the break-up stage were taken into 
account [23]. The transition energy E™* for the begin- 
ning of restoration of the properties of isolated nuclei was 
chosen to be 1 MeV/nucleon. Such value has, however, 
to be considered as a lower limit [23|. 

The isoscaling and the isobaric yield ratio methods 
were applied to the SMM-calculated primary fragment 
yields including all the isotopes for each element in the 
calculation. The symmetry energy coefficient to temper- 
ature ratio, C sym /T, was extracted for each A from the 
isoscaling parameter a and the isobaric yield ratios ac- 
cording to Eq. [T] and Eq. [21 respectively. To allow a 
comparison of the two methods, the results C' sym /T may 
be plotted as a function of A for the isobaric yield ra- 
tio method and versus A = 2Z for isoscaling. For the 
isobaric yield ratio method fragments produced in the 
78 Kr source were considered. Similar results were ob- 
tained from the 86 Kr source. The symmetry energy coef- 
ficient extracted from the data analysis with both meth- 
ods will be referred to as C°^,, to distinguish it from 
the input symmetry energy coefficient, C™ TO . The micro- 
canonical temperature predicted by SMM was used as the 
source temperature, T. In Fig. [1] the weighted average 
of the C°ym values obtained for each A from the isoscal- 
ing (panel a) and isobaric yield ratio method (panel b) 
for different E* are plotted as a function of the input 
Cly m value. The expected values, C™4 = C™ m , are 
also plotted for reference (full red line). The error bars 
reflect the statistical uncertainties in the fragment yield 
and in the temperature. The C°™^ values are monoton- 
ically increasing with C™ m , independent of the method 
used, but the calculated values are systematically higher 
than the input values. The discrepancy could be as- 
cribed to the value of the temperature used to deter- 
mine C°"* n . Indeed, the microcanonical temperature is 



o 10 



(a) Isoscaling 
~-r in -r OLJI 

" OE*=3MeV/nuc y 




(b) Isobw 


-AE*=5MeV/nuc 






_ 0E*=7MeV/nuc j^*^^ 













10 20 10 20 



CSJMeV) C^JMeV) 

Figure 1. (Color online) Correlation between Cly m and C°y^ 
values obtained from the analysis of primary fragment yields. 
C°ym is obtained by the isoscaling (a) and the isobaric yield 
ratio (b) methods for all the E* and Cly m combinations. 

different from both the "kinetic" [24| and the "chemical" 
[25l | temperatures that can be extracted from particle ki- 
netic energies and isotopic yield ratios, respectively, and 
that are typically used in similar analysis [y, [HI, [26[ . Es- 
sentially, no E* dependence of C"^ values is observed. 
A small deviation of the order of 10% is present only 
for the isoscaling-extractcd C°f m at Cf ym = 25 MeV for 
E* = 7MeV/nuclcon. Comparing panels a and 6, we 
observe that the values extracted with the two meth- 
ods applied to primary fragments are similar, though a 
slightly larger deviation from the expected values is ob- 
tained by the isobaric yield ratio analysis rather than 
the isoscaling. Moreover, the observed correlations sug- 
gest that conclusions on the symmetry energy coefficient 
value, with a systematic uncertainty on average of ~ 20% 
and 40%, could be drawn by analysing primary fragments 
with isoscaling and isobaric yield ratio methods, respec- 
tively. 

Several works aimed to reconstruct primary fragments 
from measured quantities @, Nevertheless, in the 

majority of the experiments such information is not avail- 
able. Therefore we analyse the impact of the secondary 
decay on the C™ m to *. correlations that can be ob- 
tained with the two methods. The same analysis pro- 
cedures applied to primary fragment yields were applied 
to secondary fragments to build the C™ m vs, C°yf n cor- 
relation presented in Fig. O Values obtained by the 
isoscaling and the isobaric yield ratio methods are plot- 
ted in the left and right panels, respectively. Note that 
the y-scalc is different in the two columns. C™4 values 
were extracted including in the calculation all fragments 
(referred to as inclusive - panels a, d), 4 < A < 12 (panels 
b, e) and 16 < A < 24 (panels c, /) fragments. We recall 
that A = 2Z for the isoscaling. We now focus on the 
values obtained with the isoscaling method (panels a-c). 
Largely, (7°™*. increases monotonically with C™ m . This 
is true both for inclusive (panel a) and light fragment 
(panel b) isoscaling-extracted values, independent of the 
source E* . The exception is the E* = 7 MeV / nucleon in- 
clusive case, in which a C°y^ of about 10 MeV is obtained 
for C"y m = 25 MeV. This indicates that high source E* 
might be an issue when extracting including heavy 

fragments in the analysis. A reasonable agreement within 
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Figure 3. C 3y m/T values determined by the isoscaling (panel 
a) and the isobaric yield ratio method (panel b), from pri- 
mary, secondary and filtered secondary yield distributions. 
The values obtained from GEMINI deexcitation and those 
obtained for Ar and Hf sources are also plotted (see text). 
C"y m = 25MeV and E* = 5MeV/nucleon were used in the 
calculation. 



Figure 2. (Color online) Correlation between Cly m and 
C°ym values obtained from the analysis of secondary fragment 
yields. C°y^ is obtained by the isoscaling (a-c) and the iso- 
baric yield ratio (d-f) methods for all the E* and Cly m combi- 
nations. In the analysis all fragment masses (a,d), 4 < A < 12 
(b, e) and 16 < A < 24 (c, f) were considered. For isoscaling, 
A = 2Z. 



30% - 40% of the value of C™ m and C°f m is observed for 
C™ m < 14 MeV. Above such value, the increase of 
with C l S y m is reduced for the inclusive and 4 < A < 12 
cases, while a saturation is observed when applying the 
isoscaling to heavy fragments (panel c). This is proba- 
bly due to the decay toward stability of exotic isotopes. 
Therefore, within the SMM framework, experimental val- 
ues of or close to 14 MeV, recently reported in [f| |28l - 
[30j . should be interpreted as the lower limit for the real 
Cly m value, especially if heavy fragments are included 
in the analysis. The stronger dependence between C™ m 
and when only light fragments are considered in- 

dicates that a better understanding of heavier fragment 
production should be sought. Reconstruction of primary 
fragments could provide a better understanding. Cor- 
relation techniques may provide a powerful tool for pri- 
mary fragment reconstruction 0, \2l\ . A small depen- 
dence on E* is observed for any mass range (except for 
C™ m = 25 MeV). This indicates that the different effects 
of secondary decay on fragments with different excitation 
energies mainly cancel out when taking the ratio of the 
same isotope produced in two similar sources 

A monotonic increase of C™4 with C™ m is observed 
for the isobaric yield ratio method (panels d-f). This is 



true independent of the mass range. A strong disagree- 
ment of a factor of up to 4 — 6 is found between the 
input and the output values with an inclusive analysis 
(panel d) and when limiting the analysis to heavy frag- 
ments (panel /). A better agreement between C™ TO and 
C°ym i s obtained for light fragments (panel e) . Although 
the output values are still higher than the expected val- 
ues, they present a behavior similar to the one observed 
for primary fragments and appear to be only shifted by 
a constant offset. This suggests that, with the aid of a 
simulation to estimate the shift, the input values could 
be retrieved. In addition, a source E* dependence is ob- 
served for the C™4 values. This indicates that, when 
analysing experimental data, a source excitation energy 
selection should be applied to the data to extract mean- 
ingful information. It is clear that, within this picture, 
the Cgym values extracted with this method are strongly 
affected by secondary deexcitation effects and it is more 
difficult to restore the correct C™ m for hot fragments 
without relying on a simulation. 

To investigate in more detail the effect of the secondary 
deexcitation on the C™^, values extracted with the two 
methods, we analyse the dependence of C S yl n jT on the 
fragment mass A. In addition to results on the fragmen- 
tation of 78 ' 86 Kr sources, fragmentation of 39 ' 43 Ar and 
I56,i72jjf sources w ith E* — 5MeV/nucleon were also 
simulated by the SMM model. The four sources have the 
same N/Z of 78 Kr and 86 Kr, but charges half and double 
the Kr charge, respectively, and were used to investigate 
possible finite size effects. 

In Fig. [31 C°y* n /T values obtained for both primary 
(triangles) and secondary (stars, squares and circles) 
fragments from isoscaling (panel a) and the isobaric yield 
ratio (panel b) methods are plotted for 6™" m = 25 MeV 
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and E* = 5 McV/nuclcon for the Kr sources. Values 
of C°y^/T obtained for secondary fragments for 39 > 43 Ar 
(crosses) and 156 > 172 Hf (X) sources are also plotted and 
will be discussed later. The secondary deexcitation was 
performed, for comparison, both by SMM (stars and cir- 
cles) and GEMINI (squares) [3l| applied to primary frag- 
ments generated by SMM. Indeed GEMINI includes in- 
formation on the nuclear structure and the energy levels 
for higher mass nuclei, in contrast to a more schematic 
description included in SMM 0. The secondary yields 
filtered for the experimental acceptance, thresholds and 
resolution are labelled as "+ Filter" (open symbols) and 
are also plotted in the figure. 

We now focus on the C°"^/T values obtained for pri- 
mary fragments (triangles) with both methods. As al- 
ready mentioned, the values obtained with the two meth- 
ods are in good agreement. Moreover the C°y* n /T val- 
ues show a slightly increasing trend as a function of 
A, more pronounced for higher C™ m . This dependence 
of Cgyln/T on A could be attributed, according to our 
analysis, neither to Coulomb contribution, nor to differ- 
ent average temperature at which different masses are 
produced. The C™^/T values extracted by isoscaling 
from secondary fragment yields (open symbols and stars - 
panel a) are lower than the values extracted from primary 
fragment yields (triangles). Moreover a constant behav- 
ior of C° y l !; n /T as a function of the fragment mass is ob- 
served, as opposed to the slightly increasing behavior ob- 
served for primary fragments. The obtained values show 
a gap between A = 14 and A = 16, not present in the 
value obtained by GEMINI, which is due to the change 
of the deexcitation mechanism implemented in SMM for 
A = 16 0. The values obtained by SMM and GEMINI 
are in good agreement for A < 16, within statistical un- 
certainty. A good agreement of the C^y^/T extracted 
from filtered (open symbols) and unfiltered (stars) sec- 
ondary yield distributions with the isoscaling method is 
observed, which indicate that the detection efficiency of 
the apparatus does not affect the values obtained with 
this method. Finally, the C°y^/T values obtained for 
the Ar and Hf sources are in good agreement with those 
obtained for the Kr sources. This suggests that possible 
size effects, if present, also do not affect the isoscaling- 
extracted values. 

We now focus on the values extracted with the iso- 
baric yield ratio method (Fig. [3] panel b). As opposed to 
the behavior observed for isoscaling-extracted C°^/T 
values, the C°y* n /T values determined from the yields 
of secondary fragments are higher than those obtained 
from primary fragment yields and present an increasing 
trend as a function of A. Within the SMM model, such 
a trend can be ascribed to finite size effects. Indeed a 
much stronger increase and higher values are observed 
for the C°yf n /T obtained for the Ar source (crosses), 
while a more constant behavior and values lower than 
those obtained from the Kr source are found for the Hf 
source (X) . This is in agreement with the observed finite 
size effects affecting the yield ratio of isobars reported 
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Figure 4. (Color online) Experimental (full symbols) and 
SMM (open symbols) C sym /T values obtained with isoscaling 
(circles) and isobaric yield ratio (triangles) methods. For the 



SMM calculation, CI 



14MeV was used. The excitation 



energy was 5MeV/nucleon (see text). 



in [32j . within canonical and grand-canonical models. 
Moreover, the values obtained by GEMINI arc not in 
agreement with those obtained by SMM and are rather 
scattered, as opposed to the increasing trend observed 
for the SMM-calculated values, especially for A > 14. 
Therefore, while the effects of the secondary deexcitation 
are mitigated when taking the ratio of the yields of the 
same isotope produced by similar sources, as suggested 
m Ref. [10(|, this may not be the case when taking the 
ratio of isobars produced by the same source, as done 
for the isobaric yield ratio method. In addition, the ex- 
perimental filter response causes an increase in C^y^/T 
for 13 < A < 29, as opposed to the good agreement 
observed for the isoscaling method. A limitation of the 
isobaric yield ratio method as compared to isoscaling is 
the more severe influence of the background, of the rel- 
ative contamination from neighbouring isotopes for the 
rarest isotopes and of the identification efficiency of the 
experimental setup. Indeed these effects cancel out in the 
isoscaling method, while this is not true for the isobaric 
yield ratio method, where the ratio of the yields of differ- 
ent isotopes is considered. Therefore, the isobaric yield 
ratio technique is more sensitive than the isoscaling to 
the experimental limitations; the C°™4 values extracted 
for large masses are particularly unreliable. 

Finally, we compare the SMM predictions to experi- 
mental data to infer, within the SMM framework, a C*™ m 
value. Recently isotopic yields of fragments produced in 
quasi-projectile multifragmentation of 78: 86 Kr+ 58,64 Ni 
at 35 MeV/nucleon were measured with the 47r NIMROD- 
ISiS array [3, Good isotopic resolution was ob- 

tained for fragments up to Z = 17. Information on the 
neutron multiplicity was provided by the TAMU Neu- 
tron Ball [33| . in which the charged particle array was 
housed. The mass, charge and excitation energy of the 
fragmenting quasi-projectile source were reconstructed 
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on an event-by-event basis, allowing a selection of a well 
defined source. The charge of the source was constrained 
to be Z = 30 — 34. Details on the experimental setup and 
on the source reconstruction can be found in Ref. [2l| . 

It has been reported in the literature [2l|, HH that 
better isoscaling can be obtained when selecting two 
sources with well-defined isotopic compositions rather 
than performing a system-to-system isoscaling. The 
quasi-projcctiles were therefore identified by their ex- 
citation energy per nucleon and their relative neutron 
excess m s = N ~, , where TV, Z and A are, respec- 
tively, the number of neutrons, protons and total nu- 
cleons in the reconstructed source. The excitation en- 
ergy of the quasi-projectile source was restricted to be 
between 3.5 MeV/nucleon and 5MeV/nucleon (E* = 
4.6 MeV/nucleon). Two bins in m s were selected, which 
had average values of mjj" = 0.097 and = 0.183 
(m^" — mjj" = 0.086, to be compared to 0.0859 obtained 
for the 78 Kr and 86 Kr sources, input of SMM). Isoscaling 
was performed between the two sources, obtaining a value 
of A = 0.0370 (see Eq. consistent within 2% with the 
A value used to analyze the SMM data (A = 0.0378). 
The isobaric yield ratio values for C sym /T were extracted 
for a source whose average m s was equal to the difference 
in m s of the sources used for the isoscaling. This allows 
us to compare the results from the two different methods. 

In Fig. 2] the experimental C sym /T (full symbols) and 
the SMM C° ym /T prediction for secondary fragments 
(open symbols) obtained by the isoscaling (circles) and 
the isobaric yield ratio (triangles) methods arc plotted 
as a function of the fragment mass A. The values pre- 
dicted by SMM were obtained for C™ m = 14McV and 
E* = 5 MeV/nucleon, which gives the best agreement to 
the experimental mass and charge distributions. We re- 
call that a value of C l sym ~ 14 MeV is in agreement with 
previous experimental results [1, [H, [2^] , and should be 
considered as a lower limit for the real C™ m . The pre- 
dicted (open circles) and experimental (full circles) values 
obtained by isoscaling method are in good agreement, es- 
pecially for masses A > 16. We observe that SMM repro- 
duces both the experimental trend and the values. On 
the contrary, the predicted (open triangles) and experi- 
mental (full triangles) values obtained by isobaric yield 
ratio method are not in agreement. This is true espe- 
cially for masses A > 16, where SMM predicts C°^ m /T 
values higher by a factor of about 1.5 — 2 than the ex- 
perimental values. Moreover the experimental data are 
approximately constant as a function of A for masses 
A < 20. The data for mass A = 21 and A = 23 have 
a larger uncertainty due to the worsening of the isotopic 
resolution as the fragment mass increases. Moreover, the 
already mentioned sensitivity of the isobaric yield ratio 



method to identification efficiency effects makes C° ym /T 
values extracted for higher masses (A > 20) less reliable. 
The increasing trend predicted by SMM is not present in 
the data, within our isotopic resolution. This might sug- 
gest that a more accurate description of the deexcitation 
stage should be implemented in the model. 

To conclude, the effect of the secondary decay as 
predicted by SMM on the symmetry energy coefficient, 
C° syml extracted by the recently proposed isobaric yield 
ratio method was investigated and compared to the sym- 
metry energy coefficient obtained by isoscaling. Mono- 
tonic increase of C° ym with C l sym independent of E* 
were observed for C° ym values extracted from primary 
fragment yields, suggesting that conclusions on C sym 
can be deduced from the yields of primary fragments, 
with both methods. A reasonable agreement of C™ m 
and C™^ was found for secondary fragment yields anal- 
ysed with isoscaling method, in particular when limiting 
the analysis to lighter fragments. Experimental values of 
C sym > 14 MeV should be considered as an upper limit, 
due to the saturation of the C° ym values. The isobaric 
yield ratio method retains a strong dependence of C° ym 
on C™ m , but the obtained values arc higher than the 
expected values by a factor of « 2 in the best case. 

The analysis of C° ym /T values extracted with the two 
methods showed that the secondary decay affects differ- 
ent observables to a different degree. In particular, the 
isobaric yield ratio method is more vulnerable to sec- 
ondary decay and detection and identification efficiency 
effects than the isoscaling method. This suggest that ob- 
servables built from the ratio of the yields of the same 
fragment should be preferred to observables dependent 
on the absolute yields of different fragments, according to 
SMM predictions. The discrepancy of SMM and GEM- 
INI results indicates that the impact of the secondary de- 
excitation should be further investigated, possibly with 
a more refined description of the secondary decay. 

Finally a comparison of SMM predictions to experi- 
mental data was presented. The SMM predictions repro- 
duce well the C S y m /T values obtained with the isoscaling 
method for C™ m = 14 MeV, in agreement with previous 
experimental results, which should be considered as a 
lower limit for the real C sym . 
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